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ABSTRACT 


In the interest of developing more thorough methods 
of analyzing ship structures, an entire bottom section 
extending from side to side and from one bulkhead to 
another is analyzed as an orthotropic shell of variable 
curvature. The longitudinally framed destroyer type ship 
is studied under hydrostatic loading and longitudinal 
bending stress. Specific consideration is given to the 
orthotropy of the structure and its curvature. 


As a result of the combined analysis of the shell 
and its associated stiffener system, partial differential 
equations are developed which define ‘the behavior of 
the stiffened shell. The inherent complexity of the 
developed equations precludes an analytic solution. An 
approximate solution could be obtained, however, by the 
use of difference equations and digital machine computation. 


An application of membrane theory to an equivalent 
orthotropic shell provides a means of determining direct 
stresses and extensional rigidities. A comparative 
study of two existing destroyer types using developed 
theory reveals that membrane theory with corrective 
iterations is of value in preliminary analysis. the 
sensitivity of the variation of hull curvature from the 
analytic curvature is also demonstrated. The use of 
the midship section coefficient as a design parameter 
for the structural designer should be encouraged and 
exploited. If the midship section area is kept constant 
the curvature can be varied to the advantage of the 
structure with no effect on speed. The relatively high 
bending stiffness in the girthwise direction of current 
designs merits a reevaluation of transverse strength 


requirements. 


Structural design procedures employing orthotropic 
rigidity coefficients in combination with either the 
solution to the above differential equations or a modi- 
fied membrane theory, have been outlined by the author. 
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NOTATION 


= transverse web frame spacing (in.) 
a! longitudinal stringer spacing (in.) 


- oO (  )! 


Y 


= external load (lb/ft) 

+ hydrostatic load (1b/ft°) 

= radius of curvature (ft) 

- shell thickness (in.) 

= equivalent shell thickness (in.) 

- axial displacement 

= girthwise displacement 

= radial displacement 

- longitudinal coordinate (ft.) 

- distance from shell middle surface 

= cross sectional area of a stiffening member (in. 
> beam of ship at design waterline (ft.) 


- dist. of centroid of stiffening member from 
shell middle surface (in.) 


- extensional rigidity (lb/in) 


= 





<= 


< 
Pee ee ke 


modulus of elasticity (psi) 
shear modulus (psi) 
draft of ship at design waterline (ft.) 


moment of a ee of stiffening member about 
its N. A. (in) 


polar moment of inertia (int) 

bending rigidity coefficient (lb - in) 
length of shell between bulkheads (ft.) 
length of shin (ft, } 

moment resultant (lb-ft/ft) 
forcomeccuitent (lb/ft) 

transverse shear force resultant (1b/ft) 


rigidity moment coefficient (lb-in) 


transverse vertical edge loading on shell (1b/ft) 


longitudinal vertical edge loading on shell (1b/ft) 


shear strain 

linear strain 

Poissons ratio 

stress intensity (psi) 

elastic yield stress (psi) 

longitudinal hull girder bending stress (psi) 
shear stress (psi) 


girthwise or circumferential coordinate (deg) 


Vi 





I. INTRODUCTION 


This analysis can be considered as a contribution to 
the continuing effort being made to produce more rational 
and less empirical ship structural design procedures. ‘The 
analysis of a structure whose loading and response character- 
istics are not readily defineable does indeed present a 
Challange. There exists, however, opportunities for improve- 
ment in the static, elastic analysis of ship structures. 
Such improvements could contribute towards a better determina- 
tion of adequate but not excessive scantlings arranged in the 
most efficient manner. Scantlings are normally selected on the 
basis of a series of approximate analyses. While the effect of 
continuity is often considered, the effects of adjoining members 
and shell curvature are neglected. Checks are made on the 
adequacy of stiffeners and stiffener-plating combinations on 
a panel basis. This artificial method of analysis, while 
adequate and sufficiently approximate in the past, cannot be 
considered as satisfactory for the future. The knowledge of 


ship seaway loadings is increasing. Advances are being made 


in dynamic loading theory. If these accomplishments are eventually 


to be used to greatest advantage in design, they must be used 
in conjunction with a realistic static analysis. 

The structural analysis of a plate or shell stiffener system 
can be approached in general by either of two basically diff- 


erent view points. ‘The first is that of a grillage beam analysis 


aye 





wherein the plating is taken to act as additional flanges to 
the stiffeners in an orthogonal network, Essentially the 
grillage beam concept as originally developed by Vedeler (1) 
relys on equating deflections at stiffener intersections and 
moment distribution. The solutions may be obtained by ortho- 
dox solutions to differential equations, trigonometric series, 
energy methods or relaxation and numerical methods suitable 

for digital machine computation. Variations in the grillage 
beam technique and its solutions have also been developed (2). 
One need only to consult contemporary technical journals re- 
lating to aeronautical and civil engineering to find a prolific 
assortment of techniques and methods for grillage beam solutions 
most of which lead ultimately to machine computation (e.g. 


stiffness matrices). 


The second basic approach is to utilize plate or shell 
theory for the given structure in such a manner that the 
stiffener system effect is incorporated into the medium. This 
is accomplished by substituting for the actual shell with 
stiffeners (frames and longitudinals) an equivalent shell with- 
out stiffeners which exhibits different elastic properties in 
the longitudinal and circumferential directions (orthotropic 
behavior). The manner in which the equivalent orthotropic 
sell is @etermin@d will be discussed shortly. Let if st@iigec 
at this point to state that the result of such a substitution 


will be to obtain a shell structure whose “general behavior 


% Numbers in brackets after an author's name indicate 
references listed in Appendix E of this paper. 
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under loading is equivalent to that of the stiffened shell. 

The orthotropic behavior ofthe equivalent shell may be. called 
constructional since it does not result from a property of the 
material but from the manner in which it is distributed. It is 
interesting to note at this point that this orthotropic approach 
is diametrically opposite to the grillage beam method (i.e.,the 
grillage beam assumes plating to act as additional flange 
material for stiffeners whereas in the orthotropic shell, the 
stiffener influence is distributed to the shell or plating). 

it is felt that the orthotropic shell concept has sufficient 
merit to be the general’ approach used in this analysis. Physi- 


cally it appears to be more meaningful in the response and 


behavior of the stiffened structure. While the task of analysis 
of such a structure is not simple, it is at least amenable to a 
mathematical investigation. Orthotropic shell theory consistent 
with the structure being studied will be utilized. 

Additional methods of shell or plate stiffener systems 
include the concept that the stiffener system acts as an elastic 
support for the loading and as such may be treated by an applica- 
tion of the beam on an elastic foundation method (3). Another 
interesting technique is to disseminate the effect of one set 
of stiffeners (frames) into the plating while other (longitu- 
dinals) set is subjected to a harmonic analysis producing influence 
coefficients for deflection or bending moment for each 


longitudinal (lh). 


This analysis will investigate the structural behavior of 


the bottom section of a ship from side to side and from bulkhead 


= 3- 





to bulkhead, The ship type considered will be a longitudinally 
framed naval warship in the destroyer catagory. smmphasis will 
be given to shell curvature, continuity of the shell between 
bulkheads, and the stiffener characteristics. Normal hydro- 
static loading combined with the overall longitudinal bending 
stress will be considered. Shade (5) developed a bending, 
theory for ships bottom structure utilizing orthotropic plate 
theory. However, a flat bottom was assumed and no considera- 


tion was given to hull curvature. 
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II. ANALYSIS 


Preliminary 


The portion of the ships structure to be analyzed is 


shown in Figure ec. 





a es 2 Identification of Section Analyzed 





The boundary conditions for this section are discussed 
in Appendix © . In brief they are approximated as follows: 


The degree of fixity can be assumed to be 100% clamped at the 
junction of the shell and bulkhead and simply supported at the 


Junction of the curved shell to the vertical or flat portion 
of the ships side. While these assumed conditions are not in 
reality present, a determination of the elastic degree of 


fixity is outside the realm of this analysis. The postulated 


conditions will be sufficiently accurate for approximation 


purposes. 


The inclusion of the huil curvature in the analysis can be 


By 





effected by the selection of an analytical expression for .the 
radius of curvation (r), as a function of girth angle (Q). 

The selection of the expression for the radius of curvature 
must be uncomplicated enough so that it will not impede the 
mathematical operations involved in the analysis. At the 

same time it must sufficiently describe the hull form for 

any given B/H ratio. A consideration of pertinent hull forms 
and analytical expressions is made in Appendix A; The result- 
ing analytical expression satisfying the above requirements 
has been found to be r = 2Hcos%; the radius of curvature which 
describes a cycloid of depth H and span 7 H. One of the sig- 
nificant contentions of this thesis is that a ship type such 

as a destroyer has a hull form which cannot be treated with 
flat bottom or flat plate assumptions. Of course if a restric- 
tion is made to only a small panel of the structure, a flat 
surface is more closely approached. This type of panel analysis, 


however, is a rather limited view of the structures behavior 


and disallows for the continuity of deformation. Moreover, not 


only is significant curvature present but it is variable around 
the circumference of the hull. An assumption of constant cur- 
vature or a circular cylindrical shell would greatly facilitate 
the mathematical analysis but would not yield results which 

It is felt that the 


are representative of the structure. 


variable curvature selected represents sufficiently the type 


of hull forms being considered (see Figs.19-23). 








Isotropic Shell Analysis 

To begin the analysis the elastic laws for the variable 
curvature isotropic shell (ships hull without stiffener system) 
will be developed. The shell element and associated stress 
system are indicated in Fig. 3. 


Equilibrium of forces provide the following relations: 
Kodir: 


© (Ny) dxrdg 2 (Nox) dddx + p,rdgdx = 0 
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Kquilibrium of moments give — 


about Py axis: 


- O(Mg) dxrdg 2 OC )dGdx + Qgdxrd® = 0 
Ox xg? "4 


OMg 1OM - = 0 
Ox — Ts ¥ (4) 


about Pg axis: 


O (M, axrag 4 O (Myx )afax - Qxdxrdp = 0 
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About Pe axis: 


“Ny grdpdx + Ngxdxrdg - 2 My ,dxad = 0 


or 


NO 


(6) 


Ng - "Gxt = “px = O 


Elimination of Q, and Qg by use of eq. (lL) and -(5) gives the 


‘following equations of equilibrium. 


oN, 19 


x ay te Be = (7a) 
1 ON og , : og ~ a Ny , PG = 0 (7b) 
ids Og Ox 6 Ox re 0 of 
1 OM 2 DM BE My OfL Wf -Pe-e Ue) 
r Ox2 rv oxag r QGdx O\r dg 
Nyg - Ngx+2l Mgx = 0 (74) 

r 
At this point there are equations and 8 unknowns. A study 


of the shell deformation will provide additional relationships. 


In the analysis of deformations the following assumptions 


are made: 


(1) All points lying on one normal to the middle sur- = 


face before deformation do the same after deformation (i.e., 


neg= 


Ae 





lect deformation due to Q,, 2p transverse shears). 

(2) For all kinematic relations, the dist. Z of a 
peint from the middle surface may be considered as unaffected 
by the deformation of theshell, but that for all considerations 
of the stress system, the stress, in the Z direction may be 
considered negligible compared to ES and Vg. 

(3) All displacements are small, i.e., negligible 
compared with radii of curvature of the middle surface and 
their first derivatives are negligible compared to unity 


(linear assumption). 





____ Section Pecoust. 





Peet Displacement of Point A 


U, V, W refer to displacements in middle surface of shell 
in directions x, 0, and z respectively. 


Una, Va, Wa refer to displacements of an arbitrary point A at 
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Sie 





a point Z from the middle surface. 


ainexoyee ie Ee 


U,= U-2Zo0W 
ox 





Wa = W (assumption 2) 


\! 


The strains at point A can be expressed by the following 


relations from Flugge (6). 


9 Us 
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The strains may now be expressed in terms of middle 


surface displacement by using point A displacement relations. 


Of vit+z)-2 dw7,2_[ du-2 dw 
C0 2{% , r og ¥ pe a TEN 





A simplification can be made at this point by considering 
the shell thickness to be very. small in relation to the. radius 


at any point. In otherwords Z is negligible compared to 1, 
if 


(l4+2/r) =1 
rae 
r+Z2 " Pr 





The strain relations now become 
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Carrying out the indicated operations and remembering that 
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Mise tune Groen Of p (rv = 2Hcos%), we obtain the following 


€ 
relations when the notations O _( )' and o _( ) are 


Ox Od 
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(EO ns v - z (wtangtw) w (8b) 
p L re a 
Ww! (8c) 


<j 
Kh Ice 

t 

NO 
ry po 


xP 


“ise 





The stress resultants of Fig. 3 are by definition 


t7 2 £/2 
N, = Se. (1 + 2/r) dz ’ Ny = itr dz 
-t/2 =-t/2 
t/e2 Be 
ag = Tag + z/r) dz ; Nox = cae dz 
=t/2 -t/2 
| (9) 
t/2 t/2 
i = (,(a+ z/r) z2dz ; Mg = - {  g2dz 
we -t/2 
BLE t/2 
Mage [Ugh tafe) cae, Moe) T paste 
=t/2 -t/2 


The terms (1+ z/r) goes to 1 however by the previous 


assumption of negligible z/r ratio. 


These stress resultants can now be found with the aid 


of the strain relations (¥ ). 
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since 


_ E_ (6 + bx) 
we g+v 


Tx¢ se Y xg 


2(1+v) 


After substitution of (|g) and ( 8) in ( 9 ) the 


following stress resultant relations are obtained where 


D. 25. is the extensional rigidity and 
—— 
Ke _ is the flexural rigiditveors>cndine 
12(1-.*~) 
Stiffness. 
Ny mg 0 logo, ALS) 5% + VD w) 
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Ng See ee > rut) (1D) 
r 


Nig = Ngx = D (L-¥) (Q + rv!) 
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These relations represent the stress resultants of 


an isotropic cylindrical shell with variable curvature 





r= 2Hcosf . (Note that the shell thickness has not been 


fixed yet but can still be a function of p). These relations 


exist in the unstiffened ships hull. 


We consider now the gridwork of frames and longitudinal 





stringers only, in Figure 5 





Sectiow K= const, ¢ 


Fig. 5 Element of a Gridwork Shell 
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N. = EAg € gy where €& g = strain at centroid 


by by 


Using the strain relations of (eqe8)with z = Cg 


Ng _ Bay eee) 
bj r re Be 
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Mg = My 2 
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by r 
note that MA x Jg ee 
Mx Jx by 


Summarizing the gridwork stress resultants of a variable 


curvature cylindrical shell form (r = 2HcosY) . 
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(12 con't) 


NOTE: Nx¢ and Nax could also be determined but it will be 


seen shortly that they will not be necessary. 


These relations exist in the ships hull stiffners. 


Superposing the stress resultants of the isotropic shell 


With those of the gridwork gives final stress resultants of 


the stiffened variable curvature shell (or the stiffened ships 


hull). 
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Where the orthotropic rigidity coeffecients are defined as 
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Note that the No and Ng, represent shear carried by 
the shell alone. The gridwork can only carry this shear 
by a bending of its stiffeners, a bending which becomes 
impossible when the shell is joined to the stiffeners. 
Introducing, now, these new stress resultants for 
the stiffened ships hull, into the equations of equilibrium 


(eqs. 7), one obtains the following differential equations. 


ey ye ar 2 | Oe) 
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After carrying out the indicated operations, rearranging 


terms and using partial operator notation as follows, 





qd, = ey ag )! 3 do. O_ . | ) 
"Ox Og 
and dy” = ye 9 do” = o> 
Se Sg 


the differential equations assume the forn, 
2 2 ( 15a) 
\e?4f + Cord -, + Cortangda | uy | Cr dy do | V 


2 Re 2 
Pye dy - Our7dy? + Cgayag” ¥ 2C, tand aya |w ~ -r-p 


2 
e574 | ui + | ox", ar Cardo* 4 Conta? + Cor tang a, 
| 2 
= Cer tang ds V + | -Cod>, = Cgrdo- : Cyor“ajd, + C7r@a, 
(15b) ji 


2 
-3Cgr tan<p d, - Cj,d9 = 30), tan @d., - Cer tang oe 


2 
-2C9 tang Le - 03, tang? ado* + C7 r tan + Cer tang] w= -r4pg 


ee 





|-esxita + 4 O27 | De | ~ Cardo? - 3Cgr tang 4,° 
ae Cord Tb tan“P do - Cgrde | V ce 
| Saray + 0, gubaylas? + Godt 4 C, 4rtang a “a, 
+ Cy) tan? do? - 2¢grd2° + 8Cotan“p do~ + 20,a,° (15¢) 
= 26,145 + 1054 tan“ de - 4Cgrtang do + 90, tan? 4, 
+ 15¢,, tan7¢ a, 4 Core - 3Cgr tan“ - Cpr | w= rtp, 


The constant part of the coefficients are identified as 


Cz = Dy+Dxg Cio = Kyt+ Kg 
c= 
ieee) wi 
D 
a Cla = K 
6 = Kg, “13 = Kyg + Bg. + Ky 
2 K 


Equations (15, a, b, c) comprise a set of 3 partial differential 


equations in 3 unknowns, (u,v,w). Their solution as a function 


2 





of radius and rigidity coefficients will completely define 


the structural behavior of a stiffened cylindrical shell 

(of cycloid cross section) for given loading and boundary 
conditions. The rigidity coefficients were made constant 

as a necessary simplification. It will be seen later, how- 
ever, that many are in actual practice essentially constant 
over considerable portions of the shell section. Additional 
comments concerning these differential equations will be 


delayed until the Discussion of Results section. 


Membrane Analysis 


Since the solution-of the developed differential 


equations is not readily available, one must approach the 
shell problem with other and more immediately useful tools 
of analysis. It must be stated at this point, however, that 


the analytical approach by means of the differential equations 


of the problem is the only means by which a complete and 


relatively exact solution may be obtained. It is unfortunate 


that their complex nature is currently limiting their use 


by shell designers. The alternate methods employed to 

obtain solutions include the beam method, iterative paiente= 
tion techniques, membrane theory, and other unique approxima- 
tion approaches. Lundgren (7) describes rather completely 
the various methods of shell analysis currently in use. 

The applicability or the relative desireability of the 
alternate methods is determined basically by the shell dimen- 
sions. The ratio of the shell length to the shell span (or 
beam), the L/B ratio, normally will govern the gecniitiue-. The 
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application of membrane theory (i.e., the consideration of 
direct stresses only; the support of the load by tension, 
cpmpression or shear forces only) is considered good for 
short shells or shells whose L/B < 1. Lundgren (7) states 
that the membrane condition of stress is a very good 
approximation for the short shell since the elastic support 


between end sections (or bulkheads) is quite stiff and 


loading is balanced by N. - Ng » and eG stresses. The 
ASCE Manual of Engineering Practice No. 31 recommends a 
Superposition of membrane stresses and those created by 
boundary conditions. However, on the matter of boundary 
conditions, the manual also states that for a L/B < .83 
the transverse (Bulkhead) end effects damp rapidly and tke 
effects of the longitudinal edges are felt only in a small 
area. 

The previous discussion indicates that a membrane 


analysis of a ships bottom section between bulkheads is a 


reasonable approach when one realizes that the L/B range 


for the ships or ship type under consideration is as indica~ 


0.05 | outa] ones | 0.56 


= length between bulkheads ( IN & veaqion) 
beam at waterline 


ted in Table l. 













L/B 


TABLE 1 Leneth to Span Ratio of Ships Considered 
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The membrane stresses (lb/ft) in the shell element can 
be obtained by setting the shell moments equal to zero in 


the equations of equilibrium (eqs. 7). 


Ng = PL? 


od 


—— Sa 


integrating the latter two of these, 


J G04 a) tx + ae 


2 - | Gey 2 Dts) dx + f5(9) 
r Og 


I 


Ng 


= 
i 


for p, = 0O and Pg and py, independent of x. (These are correct 
conditions for ships hull loading. The longitudinal stress 
resulting from the hull girder action and determined from 

the longitudinal strength calculation will be accounted for 


shortly but will not be considered as a p, loading at this 
pommG 


- (Pg 4 1 d Xp x + f3(¢) = ~xF ig) + f,(9) 
rd? Where FO) - (Pep + 3 aNe 





Ne 2 xe dF(¢) - x df; (9) i f5(9) 
2r ag r dg 
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F(Q@) is determined by loading and shell shape 
(curvature) but f,(@) and f5(%) will be determined by 
boundary conditions at the shell ends (bulkheads). 

Consider first the case of simple support end 
(bulkhead) conditions. 


At X=+41L/2, WN, = 0 








then, Fig. 6 Shell on Simple Supports 
te aF(O) - L 4f\(0) f2(9) 
= —— = _ a aa ts 
Br ag er dg 


L° =a F(9) : L aty (9) Fat?) 
Br ag 2r dg 





which gives, 


£3(f) = 0 fo(p)= - LE aR(p) 


Br dg 





and therefore, 


Nig = - xF(O) 
N, = -1 (L° = yx®) aF(9) 
Br dg 
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Now consider shell end or bulkhead conditions for the 
shell continuous over the supports and assume 100% clamping. 
As an aid in determining the stress Ny over the supports, 
one can make use of the relationship given in ASCE Manual 
No. 31: the ratio of the stress (in a shell over continuous 
supports) at the point of support to the stress in simply 
supported shell at the midspan is equal to ratio of the 
moment in the continuous shell at the support to the moment 


in the simply supported shell at the midspan. 





or : 
and, for a strip of shell, “oO 
Ve Me M, = -1/12 whe 
To — M a ae 
N 2/3 M = 1/8 wht 
J - O 
Nxo 
Neg - 2/3 No Fig. ( Comparison of Continuous & 
Simple Supports for Shell 
and N50 was found previously to be - Le dF (@) 


Br a 


Nat support (X =+L/2) is 4 L° daF(¢) for a shell 
l2r ag 








continuous over supports. 


Fig. 8 Shell Continuous over Supports 
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Lé ar(d) Le ar(g) - L  aty(f) _ £2(8) 
“lar ag Sr age 2r sa 

L2 aF(¢) L¢ a4F(d) Lb  afi(d)  £5(g) 
— = —— 

12r ag 8r ag er ag 


which gives, 


2lir ag 
then, 
N, = - ite - 12x°) adF(¢) (16) 
2h r ag 
Ng = Py i 
where, | 
F(O)= (pg, 2 ANg) 
r d 


Equation 16 represent the general form of the membrane 


stresses in the ships bottom section of radus of curvature 


r. 
These stresses will now be developed in a more final 
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form by considering the hydrostatic loading and the 

cycloidal ie section cylinder. 

Let qq = ¥ ho where hj= an 
arbitrary design head to account 
for wave forces on freeboard surface 


eit saod Bly 
z = H/2 (\|+cos2@) from appendix (A ) 





Ma fean SS, Loading on Ships Hull 


ey ie w= internal loading of shell 
oS Y (bot?) due to machinery weights (lb/ft) . 
q = X {he +H, (+c 0s2@) | 
Pr= -(q - w cos @) 
ee - t | hy + H/2(1 + cos2q) | + wcos @ 
Pg = - w sin g 
Pp, - 0 


Ng = = | hy + H/2(1+ cos2d) | + w cos of 2H cos @ 


Ng = -¥} ata Ho cos p 4. He cos? g | + 2wH cos“¢ (17 a 
and since F(Q@) = PO 4 1 d Ng 
r qd @ 


F(@) 


\l 


Y | bo tang ie “ sinag | - 3w sing 


N, and Nig are then found by applying equation (16). 


Ng = 7% | ¥ (no tang +3/2H sina) - 30 sing} (17) 


aa 1 2. 3 had 3g 3H cos2g | (17c) 
NL ae | ¥ (no seo 4 a ) 
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Equation (17) are the membrane stresses for the ships 
bottom section. With the aid of these relations, the necessary 
shell thickness (that is, the equivalent shell eateries s for a 
hull with no stiffeners) can be found as it varies around the 
girth. From this equivalent thickness, extensional rigidity 


coefficients (Dx Dg : Dig» Dy ) may be computed. For the 


& EB & EB 
single equivalent thickness derived there results the fact 


that D, Hil ee they both equal te or 1.1lt, see appendix D. 
EE : lL =+ 
However, in practice these coefficients are very nearly equal 


as can be seen from Figure 26. 
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Fi e 10 Direct Stresses on Equivalen 
Shell Element, | 


Where FF = the hull girder longitudinal bending stress as 
determined from the longitudinal strength calculation. 
Using the maximum shear failure theory and in particular 


the Tresca version, 


2 2 
Oy _Vinax _ (& “§% ) Tx £ on 
2 2 





(19) 
y? =( - 2. 7p +94 LY. ¢ 
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Substituting equation (18) in equation (19) and 
solving for t, gives 
ho = BY: ao TN 
2A (20) 


where A= 1)h ( (ee - re 


B= 2h (Vy , = Ny ) 
c= NS = eN Ng + Ng ve ung 


It now becomes necessary to make a judgement as to the 
treatment of Aa and = 

ce will be assumed as 33,000 psi and a safety factor of 
1.65 will reduce its value to 20,000 psi. 


CF will be approximated by Tobins relation that 
SL h L. = length of ship in ft 
(F _\j/ lg Where Ls = length of ship in ft. 
1 max = 


ane Ton/in° 
If it is desired to have V7 as a function of ship depth 


(or girth), with Fi =- 9. at the keel (f = 0) and 
1 max 


OF Me max at the waterline or N.A. (@ Eh /2)) (to allow 


for the ship being in the heeled condition), 
then Oye 3 L, (5+ 3 cos2g) (21) 
> 
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The membrane analysis has now been sufficiently developed 
to make a test of its applicability to specific examples. For 
a given ship (or design), the equivalent shell thickness 
variation can be found by means of equation (20) which requires 
a determination of membrane stresses equation (17) and hull 
girder bending stress estimate equation (21). The extensional 
rigidity coefficients are then readily calculated since the 


equivalent shell has no stiffeners and D, 
Bef et) > 





The calculated values may then be compared to the actual 
values on the ship selected for study. It should be noted 
that in the calculation of the membrane stresses, the longi- 
tudinal section (X — constant) selected for stress determina- 
tion was the section at the bulkhead. This was a logical 
design decision since a close examination of equation (20) 

(or a few comparative calculations) will indicate that analysis 
of this section results in the greatest te for any given g 
value along the shell length. In addition, the value of w or 
internal loading was assumed equal to zero. Not only will 
this be the more severe loading condition, but it is also the 
one which current design (series of flat panel analyses) 
employs. This assumption will assist in making the comparison 
more valid. The retention of this quantity in the general 
theory, however, is recommended. Additional comments on 


this subject will be made in the discussion section. 
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III. RESULTS 


The analysis of the orthotropic shell by superposmelen 
of the isotropic shell and a stiffener gridwork (with all 
effects of bending and torsional moments in the shell con- 
sidered) resulted in a set of 3 partial differential 
equations in 3 unknowns (u, v, w). They are equations (15) 
on page 22. Due to their lengthy nature they will not be 
repeated here. The poraiedents in the partial differential 
equations are variable in that they contain the radius of 
curvature in powers up to the fourth. Equation (15c) is of 
the fourth order while 15a & b are third order. Equation (15a) 
could be made homogeneous by setting py = O but equations 
15b,c must remain non-homogeneous. All three partial 
differential equations are linear. 

The membrane analysis produced more tangible results 
Lal ney fom of expressions for the membrane stresses as a 
function of: girth angle (f), position along shell length (X), 


and ship characteristics such as ship length (L,), draft (H), 


compartment internal loading (w), distance between bulkheads 


\ 


(L), and a design head on freeboard plating (ho). : 


Membrane stresses in ships bottom section.* 
ee Le = 12x° ) 36 ~ Y (ho sec + 3H cone 
48H =a cas 
Ng - ¥ (2H No cos + 2H“cos30) + ew Heos“p 
Nig - my (hotang + 3/2 H sin20) - 3w sing | 
* A bottom section constructed with no stiffeners but an 
equivalent thickness (te) to produce the same extensional 


rigidity as one with stiffeners. 
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Figurell illustrates the membrane stress distribution 


along the length of the shell (a section ®= const). 


Figure \2 represents the circumferential distribution 

of membrane stresses in DE 1033 and DD 931 hulls.(a:séction x const). 
The derived expression for the shell equivalent thick- 

ness variation utilizes the membrane stresses along with 

the hull girder bending stress ((@}) and a yield or 

allowable design stress (Gy). 


: + e : 
Os es B - Bo = AC inches 


2A 
where 
A= bh (Go - 67) >» Cy» = psi 
3 Gi Ey (N, - Na) membrane stresses = lb/ft 
C = NS - 2N Ny + Neg + LN 


The extensional rigidity coefficients are then computed 


by D Dg ee 1.1lte org as (OS 


. “— 





5 1 ye 


Did te _ S20). 





D, _te | +3 3te 


en) 


Using the DE 1033 and DD 931 futiejeceese coefficients 
have been calculated from the developed expressions above and 
compared to those calculated from the actual hull structure. 


The results are plotted in Figuresl,15,16,17, 
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Figure 11 Longitudinal Distribution of 





Direct Stresses in Ships Hull 
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Figure 1 Extensional Rigidity Comparison 


for DE 1033 
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Figure 15  Extensional Rigidity Comparison 


je 


for DE 1033 
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Figure 16 Extensional Rigidity Comparison 


for DD 931 
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Figure 17 Extensional Rigidity Comparison 
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IV. DISCUSSION OF RESULTS 


The inherent complexity of the developed set of 
partial differential equations precludes a simple solution. 
An analytical solution appears to be a formidable challange. 
Perhaps the most promising course of action would be to 
apply the technique of difference equations and with the 
aid of a digital computer obtain an approximate solution. 
The complex nature of the equations should not discourage 
efforts to solve them, for in their solution lies the 
complete behavior of the stiffened shell. This last state- 
ment needs qualification in that the behavior does not 
inelude the effects of concentrated loads, or localized 
stresses. However, as noted and discussed in Appendix B 
the general behavior of the shell should be obtained first 
and then a follow up check made for localized stresses. 

The shell analysis should not proceed solely on a series 
of localized stress analyses throughout the structure, 

The value of the above equations is naturally dependent 
upon the validity of the included assumptions. Perhaps the 
most questionable one of these is the assumption of con- 
stant orthotropic rigidity constants with variation in g. 
This simply states that the structure is uniformly con- 
structed of the same size structural elements around the 
section. This, of course, would never be done since the 
loading varies with d. Figure 26 illustrates how the 


actual coefficients do vary with Q. It can be seen that 
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there are regions of the shell where many are constant 
with @. For those that vary slightly, a mean value might 
be taken. The sharply varying ones are associated with 
the stiffners in the @ direction. Perhaps a reanalysis 
of the transverse strength requirements will reduce 

these to more smoothly varying values with Q. 

Once the equations are solved, the solutions will be 
in the form of expressions for the deformation as a func- 
tion of X, 0, and the rigidity coefficients at any given 
point (with, of course, the effect of the injected boundary 
conditions being incorporated in the functional expression). 
Once these deformations are known, expressions can be 
obtained for the strains, moments, shears and stresses in 
the equivalent orthotropic shell. The rigidity coefficients 
may then be determined at any point by assignment of design 
stresses. These design stresses would, of course, be 
influenced by considerations of hull girder bending stress 
and stability of the structure. After all the rigidity 
coefficients are determined (i.e., the structural member 
sizes and spacings) they could be fed back to the original 
expression for deformation along with their variation in g. 


A recycling of this process until the coefficients do not 


change appreciably will give the final structural design. 
A follow up check for localized stresses and concentrated 


loadings will complete the process. 


The application of this theory and procedure to a photo- 
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elastic scale model and later a larger steel structural 
model would prove invaluable. A program of testing, 
modifying theory and reorienting basic approaches would 
most probably result in new and improved concepts in 
naval structural design. 

A few comments are in order eoncerning the lavecdhal 
curvature. The reasons for selecting the analtical 
curvature (or radius of curvature) used in this thesis 
are amply covered in Appendix A. However, for additional 
research in this area, NACA TM 1302, pp. 8 - 16 is 
strongly recommended. The topic of variable curvature is 
discussed and fourier po amee i are developed for the 
radius of curvature expression. Very simple expressions 
result from the analysis and they could be used to approxi- 
mate, very closely, any given hull design. It is quite 
interesting to note that within the scope of preliminary 
hull design, there exists a design variable unused by the 
structural analyst. This design parameter is the miship 
section coefficient (the ratio of the immersed area of 
the midship section to the product of beam and draft, Am). 
This coefficient defines the relative fullness at the _ 
midship section to a circumscribing rectangle. In effect 
it is a shape factor which determine the shell section 
profile. If the stress analyst had some degree of control 
over this coefficient, a shape or section profile could be 


determined which would be the most beneficial for the stress 
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in the hull. the stress magnitude and distribution in 

any shell structure is a function of its shape or curva- 
ture. ‘he encouraging fact is that indeed the analyst 
can play some part in determining this coefficient to 

the advantage of the structural requirements without 
affecting the resistance of the ship. Db. W. Taylor (8) 
demonstrated this fact in 1908 when he conducted model 
experiments to determine the influence of midship section 
shape upon the resistance of the ship. His tests were 
definitive enough for him to make the following conclusion; 
"for vessels of usual types and of speeds in knots no 
greater than twice the seaene root of length in feet, the 
naval architect may vary widely the midship section 
fullness with no material beneficial or prejudicial effect 
upon speed." It must be remembered, of course, that the 
above statement applies to the shape and coefficient of 
the midship section for a given area and not to the area 
itself. Taylor varied the midship section coefficient 
from 0.7 to 1.1 using 10 models and holding the midship 
section area constant on all tests. Thus it appears that 
the structural analyst can accept a given design and have 
considerable freedom to determine his optimun shell shape 
(for a given midship area) within the limits of not adverse- 
ly affecting the Prismatic coefficient or general 


seaworthiness. 


The use of membrane theory for the evaluation of stress 
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in the ships bottom section has considerable merit. 
Although the stiffened ship structure will carry the 
stresses according to the distribution of its members, 
membrane theory will predict what the stresses are that 
must be carried. There are two basic limitations to the 
theory. (1) The normal shears, bending moments and 
twisting moments are neglected, (2) the membrane theory 
becomes invalid at and adjacent to the boundaries. These 
limitations are not severe, however, since the quantities 
of normal shear and bending and twisting moments are 
normally of a very minor magnitude. (It will be shown 
shortly how this can be determined quantitatively). The 
boundaries can be subjected to a separate discontinuity 
analysis and do not jeopardize the evalution of conditions 
over the major portion of the shell. To prove that the 
membrane theory gives is a good representation of the 
situation, one can proceed in the following manner (1) find 
the displacements associated with the direct stresses pre- 
dicted by membrane theory; (2) use these displacements in 
the general equations which incorporate bending (equations 10); 
(3) compute the direct stresses from these bending theory 
equations and compare them to the membrane direct stresses 
to evaluate the accuracy of the membrane theory. 

Figuresll, le, represent the membrane stress distri- 
bution in the longitudinal (@=- const) and girthwise (X= const) 


directions respectively. Jt must be kept in mind that 
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these membrane stresses are invalid at or adjacent to 
any of the shell boundaries (bulkheads and longitudinal 
junctions with vertical ships side). A matching of 
reaction forces at these locations is necessary as well 
as the determination of secondary bending to satisfy 
compatibility. 

The most interesting feature of Figurell is the 
presence of longitudinal compression (- N, stress) in 
the shell midway between the bulkheads. If the ship 
actually had a flat bottom the stress would obviously 
be a tension. As the bottom curvature increases and 
becomes a shell, the longitudinal stress does indeed 
become compression. The direct stress whose magnitude 
and direction is the most sensitive to shell shape or 
curvature is the longitudinal direct stress, N,. ass eae 
not a conclusion of this thesis but rather a standard fact 
of shell structure behavior (a fact which is largely 
responsible for the widespread success and popularity 
of concrete shell structures). This is of interest to the 
naval structural engineer since it is this particular 
stress which is combined with the hull girder bending stress. 
The recognition of this stress should render ship structural 
stability considerations more realistic. 

Figure 13 represents a plot of the longitudinal stresses 
as measured on the destroyer HMS ALBUERA (9) in the hogging 


and sagging conditions. The test section (fr. 62%) was 
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located in the midship region and about halfway between 
transverse bulkheads. Considering the observed stresses 
in the hull above the waterline, a linear relationship 
is readily established. When the underwater portion of 
the hull stresses are examined it can be seen that they 
exhibit a variance from the establishedlinearity (the 
linearity established the upper portion of the hull must 
be the same as for the lower since nether moment nor 
section inertia has changed in going to the lower portion). 
This departure is obviously due to the effect of the hydro- 
static loading on the underwater shell or hull. It will 
be noted that this effect augments the negative hogging 
stress and reduces the positive sagging stress. In other 
words the hydrostatic loading on the hull produces a 
longitudinal compression (-N,). This agrees with the theoreti- 
Gad plot of Ny. in Figure 11. The longitudinal compression 
effect under hydrostatic loading would be no surprise to the 
civil engineer who works with shells. It would, however, be 
unexpected by a naval structures analyst who limits himself 
to the study of small panels and predicts only a flat plate 
tension effect due to hydrostatic loading. While it is quite 
true that each panel as such is subjected to a tension pro- 
WMcime load. this is only a minor ADORE Lon on the more 
general shell behavior which is a compression. These local- 
ized stresses can be checked in the later stages of design. 
Figure 12 illustrates the distribution of direct stresses 


around the girth (X= const) for two selected destroyer types. 
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The section portrayed is at the left transverse bulkhead of 
the basic shell (X= -L/2). This particular section was 
selected for purposes of obtaining a design shell thickness 
(and its variation with ~) since its thickness requirement 

is greater than any other intermediate section for any 

given 9 value. This was an unfortunate choice, however, as 

a forthcoming discussion will indicate. The direct stresses 
of Figure 12 were used in calculating the equivalent shell 
thickness and hence extensional rigidity constants (tables  , 
6, ). These coefficients are compared with those actually 
existing in the ships as a result of their construction. As 
was previously stated, fie membrane theory will give reason- 
ably good results everywhere except the boundaries. This 
unfortunately was the necessary region (trans. BHD) to 
omalyze sinee it produced the controlling equivalent shell 
thickness. The evaluation of the direct stresses could have 
been and should have been improved by using the procedure 
indicated on page }\8. Since the extensional rigidity coeffi- 
cients were calculated from t, which in turn was calculated 
from the direct stresses; additional corrective action on 
the direct stresses produced would result in a different and 
better correlated set of extensional coefficients. It is 
surprising that at this worst possible location for membrane 
theory, the degree of agreement between constructed and cal- 
culated is still reasonably good. A corrective cycle on the 
obtained membrane stresses would result in a general increase 


of varying degree for the rigidity coefficients at all values 
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of @ (since all Y values are on the boundary). The differ- 
ence in relative variation of calculated from constructed 
values in the two ships selected is due to several factors 
which characterize each ship. It is believed, however, 

that one of the major causes is the fact that shell curva- 
ture of DD 931 (considered similar to DD 952) resembles a 
cycloid to a lesser degree than does DE 1033. This sugg- 
ests the importance of selecting the most accurate analytical 
model for hull curvature. Other causes for non similarity 
of the ships plots are possible opposite departures from 

the assumed. conditions in the membrane equations (e.g. Fr 
stress estimate, factor of safety used and degree of fixity 
at ends). Figures 15517 > merely represent the comparison 
of the calculated and constructed ship hull SRO since 
these rigidity coefficients are just a product of a constant 
and the thickness. 

The approach used in this thesis was not expected to 
produce any radical changes in structural design. It was 
intended to achieve a more rational design attitude by treat- 
ing a section between bulkheads as a shell and then combining 
these results with the overall ship bending stress. An 
expected result was a slight reduction in structural members 
as determined by current design. With due regard to the 
previously discussed obstacles and handicaps involved, the 
plots of Figure 1)4,15,16,17, indicate an apparent trend in 


this direction. The determination of strucutral members can- 
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not be made on extentional coefficients alone. A recycling 
of the membrane solution as noted on p will give magnitude 
of moments in the equivalent shell. The individual structural 
members for the ship may then be determined by a combined 
consideration of the extentional rigidities required, the 
moments in the shell and the localized panel stress. 
Unfortunately neither one of the basic approaches used 
in this thesis (i.e., shell bending theory by differential 
equations and membrane theory) has been followed through 
sufficiently to give definitive quantitative conclusions. 
The investigation, however, has been quite profitable in 
that the ships hull has been analyzed in a manner more 
appropriate to the geometry of the structure. Such an analy-~ 
sis can and, indeed, did reveal information unobtainable 
by an otherwise less realiztic approach. Specific areas for 


additional study and analysis have been clearly indicated. 
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V. CONCLUSIONS 


1. The resulting partial differential equations of the 
orthotropic shell behavior are complex and even after 
simplification are not amenable to an analytic solution. 
2. A rational structural design procedure can be 
established by using the partial differential equations 
solution and the orthotropic rigidity coefficients (see 
page 49). 

3. The analytical expression for the hull radius of 
curvature must be uncomplicated and yet representative 
of the actual hull form being analyzed, 

kh. The midship section coefficient or the midship sec- 
tion shape is a design parameter which can be used within 
reasonable limits to benefit the structural design with 
no prejudicial effects on hull resistance or ship speed, 
5. Subject to its characteristic limitations, membrane 
theory can be used in the preliminary analysis of ship 
structures. Iterative corrective procedures can be 
used for improvement of solution and/or estimation of 
existing error, 

6. A longitudinal compression results from hydrostatic 
loading on a ship hull that possesses curvature. ‘The 


magnitude and distribution of the compression around the 
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girth of the ship is sensitive to the hull shape or 
variation of curvature. 
{. Useful extensional rigidity coefficients can be 
determined by membrane theory when iterative corrections 
are included. 
8. Current design methods produce hulls which have a 
gerthwiise bending stiffness which is five times as 
large as the longitudinal bending stiffness. 
9. The different and more appropriate approach used 
in this thesis for the analysis of ship strictures with 
hull curvature has produced; 

(a) A specific knowledge of items within this 
analysis requiring further study and refinement. 

(b) A more discerning attitude towards the 
structure being designed. 

(c) Suggested design procedures utilizing the 


methods and techniques discussed in this study. 
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VI. RECOMMENDATIONS 


1. Computer solutions to characteristic partial 
differential equations of the orthotropic shell structure 
should be devised, programmed and tested for merit in 


relation to current or other proposed design concepts. 


a. Optimization of ship structural members and their 
arrangement should be on the basis of an orthotropic 
shell and not on a series of panels or current design 


approximations. 


co A study should be conducted to determine if a 
defineable relationship exists between the longitudi- 
nal and transverse strength requirements of any given 
design. The current ratio of five to one on girthwise 
flexural stiffness to longitudinal stiffness appears 


excessive. 


he Experimental verification of all developed theory 

and techniques should proceed concurrently to aid in 
providing realistic results which are of a practical 
significance. Experimentation could commence economically 
by photoelastic stress analyses of small scale models, 


Tests on successively larger steel scale models using 
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adequate instrumentation should proceed if early results 


indicate that further testing can be of value, 


5. A short study should be made to determine an analyti- 
cal expression for the radius of curvature of any given 
design. The expression must be uncomplicated or capable 
of being simplified for the mathematical operations 

which must be made on it. (NACA TM 1302 pp. 8 - 16 is 


recommended for this study). 


6. Since the mathematical analysis of a circular 
orthotropic shell is quite complete, one possible approach 
is to break the shell into a finite number of segments 
represented by constant radii of curvature (circular), 
make the analysis and join the segments by making them 


compatible at all junctions. 


lic A strongly recommended topic for study is the 
variation of midship section coefficient or shape for a 
given midship section area and the determination of the 
optimum shape for structural design. The only major 
constraints will be the maintenance of the designed 
prismatic coefficient and desired seaworthiness character-~ 
istics. Within these constraints, however, there is 


latitude for increasing structural efficiency. 


8. A modified membrane theory should be developed for 


a5ye 





the preliminary analysis of ship structures. Since the 
solution improves with iteration, it is quite amenable 


to machine computation. 


9. A method should be devised for readily determining 
the magnitude of longitudinal compressive stress 

which is produced by hydrostatic loading on all hulls 

with curvature. This value should then be included 


in any deliberations on structural stability. 


10. A complete analysis of the stability of the 
orthotropic shell between bulkheads should be combined 
with this thesis to provide the complete ingredients for 


anew concept of ship structural design. 


ets A relieable solution to the partial differential 


equation of the orthotropic shell will require a know- 
ledge of the.shell boundary conditions. It is apparent 
that a study devoted soley to this subject would be 


of considerable value, 
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VIL. APPENDIX 
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A. HULL CURVATURE 


The curvature of the hull in the girthwise or circum- 
ferential direction is one of the most significant factors 
in this investigation. The author feels that the magnitude 
and variation of hull curvature for ships of the destroyer 
type has a significant effect upon the structural behavior 
and therefore should be given full consideration in the 
structural analysis of these hulls. The section of the 
ship considered was in the midship section region. There 
will normally be no longitudinal or fore and aft hull curva- 
ture in this middle penton of the ship. 

The argument could be made that the hull girthwise 
curvature is really negligible if a basic panel (bounded by 
frames and longitudinal stringers) is selected for analysis. 
the weakness of this Bye umienie a ereaen is that the hull's 
general response is not on a panel basis but on a much larger 
and continuous region (i.e., the portion of the hull structure 
between transverse bulkheads) wherein any flat plate or 
flat bottom assumption would be unrealistic. Additional 
comments on the relative size of hull portion analyzed will 
be reserved until Appendix B. 

The task of introducing the hull curvature of a specific 
type of ship, or of a given design within a type, is a dunn 
cult one. The analytical expression of curvature selected 
must be representative of the design being studied. It must 


also beareasonably uncomplicated function capable of being 


=OH0Q= 





subjected to the necessary mathematical operations without 
excessively expanding the complexity of the resulting 
analysis. The representation of ships hulls by ‘ olaanemaian 
expressions is not appropriate since the desired expression 
is one giving the radius of curvature as a function of the 


girth angle (@). 
de 
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Figure 18 Cycloid Profile 


The analytical expression finally selected as satisfying 
the above requirements was the equation of the cycloid 

r= ha cos g (Figure ) where the depth of the curve 
(draft) is 2a and the span (beam) is 2a. ("a" is simply 
the radius of the generating circle). This type of curve 
or section profile results in a span to depth ratio Poe B/H 
ratio) of TV . ‘The B/H ratios of the ship type considered 


is indicated in Table 2. 


a6l= 





DE 1033 DD 952 DLG 9 DLG(N) CG(N) CA 13h 


B/H 


Sh he Sr 2 alt 3.02 





TABLE 2. B/H Ratios of Ship Types Considered 

The T B/H ratio of the cycloid does not in itself 
qualify this curve as being a good representation of the 
destroyer hull. The curvature and its variation with g 
must correspond sufficiently close to the ship hull shape. 
This is the most significant factor, It must be remembered 
that coincidence of the analytical expression and the ships 
hull is not of importance (i.e., the r, © relationship should 
not be misconstrued to be a polar plot). The relative 
‘equality of the values of the radii of curvature at given 
@ values is the significant measure of the expression's 
adequacy. 

The parametric equations of the cycloid in @ are, 


y = a(2g + Sin2p) (22a) 


Cos? 
eS. + KP (22b) 


The cycloid section can then be fitted to any given design 


poeteteing a =- H/2 and plotting 


you a [oman cos7— (=) obtained from 
a a 


an elimination of g in equations 22. These plots have been 


made in figures 2Q 21, 22, 23, along with the plots of the 


curve produced by fitting an elliptical section to each actual 


An 








hull shape selected for comparison. 

A study of the comparative plots of hull shape and 
analytical expressions reveals that the ellipticel section 
possesses a greater resemblance to the curvature (or radius 
of curvature) of the actual hull than does the evel Ol amimote 
how the DE 1033 curvature is virtually coincident with the 
elliptical section). However, the analytical expression 
for the radius of curvature of the ellipse is 


acb* ; 
r—_ where a and b represent the major 


a2sin?G + p@cos“G 
and minor axis respectively. This expression for the radius 
would not be suitable for the mathematical operations in- 
volved in the shell theory of the hull. It is also of 
interest to note how the hull form of the CA 13h is not 
appropriate for this analysis since it deviates excessively 
from curvature given by either expression. CG(N) hull 
form would also be considered unsuited for the cycloid approxi- 
mation. 

The obvious disadvantage of the cycloid section is the 
unfortunate fact that the radius of curvature goes to zero 
at g elle er, This is inherent in the nature of the cycloid 
curve since at the extreme points of the span the generating 
curve changes the sign on the radius of curvature and, at 
this point, it passes through a zero point. This is not a 
major disadvantage, however, since it occurs at the region 


where the curved shell joins the flat plate portion of the 
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hull (the boundary of the shell section). This particular 
region requires a separate analysis regardless of shell 
theory employed (bending or membrane theory). In the | 
former case, rather undefineable boundary conditions exist 
along with compatibility obstacles and in the latter case 
membrane theory never does account for conditions at the 
boundaries. Therefore, it is of no conséquence that the 
radius of curvature expression given by the cycloid is 
invalid at the boundary of the shell since the shell theories 
are also in question at that region. A pertinent subject 
for investigation would be the study of the discontinuity 
stresses present at this junction, Results could be 
correlated with this thesis which purports to treat only 
the bottom section of the ship not at or adjacent to the 


boundaries. 
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B. ORTHOTROPY OF SHIPS HULL 


As a ships hull is fabricated, the elements of 
construction (shell plating, longitudinal stringers, 
transverse frames) are arranged in such a manner that the 
resulting structural elastic properties (extensional, bend- 
ing, shear and twisting rigidities) are not equal in all 
directions as they are for an isotropic structure. The 
varying of elastic properties with direction, or the aniso- 
tropy of the structure, is somewhat simplified by the 
symmetry of construction in the longitudinal and circumferen-~ 
tial directions. This special case of anisotropy where the 
elastic properties are different in mutually perpendicular 
directions is termed orthotropy and such structures are called 
orthotropic. It should be noted at this point that the basic 
material of construction is isotropic. The disposition of 
this material creates the anisotropy and therefore the 
structure can quite correctly be called constructionally 
orthotropic. 

The reaction of a structure to any given loading will 
be a function of its elastic properties. It would seem 
then that a realistic analysis could be made if the structures 
elastic properties and their change with direction are recog- 
nized. This can be accomplished by considering an "equiva-~ 


lent" structure without stringers and frames but whose elastic 
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properties and general behavior under load resembles 
sufficiently the original structure. 

The advantage of this approach is that the continuity 
of the structures behavior is retained and the resulting 
analysis is superior to one in which conservative attempts 
are made to account for the complex interactions existing 
between structural elements such as frames, longitudinals 
and shell plating. A critical question arises at this 
point - what is the characteristic mode of deformation of 
the stiffened hull? Does it deform under load to form a 
smooth half wave jenttn between bulkheads as in Figurecha 
or is all the deformation confined to the shell surface 


between stiffening members as in Figure 2hb ? 


Bey 
‘eo a 





24a. 24b. 


FPige 2 Stiffened Hull Modes of Deformation 
The answer is clearly a function of the number and size of 


the transverse frames. Lundgren (7) is satisfied to accept 
the former version of deformation (Figure 2a), whenever the 
number of transverse frames between bulkheads equals or 
exceeds four. There are panel advocates, however, who 
insist on the latter mode (Figure 2ljb) as being the predomi- 
nant one which accounts for most of the strain energy of de- 
formation regardless of the number of frames. 


For the destroyer type having approximately 4. frames 


= 7a 





between bulkheads (in the middle half of the ship), the 
author feels that the normal mode of deformation is probably 
a combination of smooth half wave length deformation with 


smaller panel deformation superposed on it as in Figurees . 





Fig. 25 Combined Modes of Stiffened Hull Deformation 


Perhaps many of those who can accept only panel deformation 
are misled by the physical appearance of ships hulls which 
have the so called “hungry horse" look or a permanent deforma- 
tion on all panels making the shell appear as though it were 
shrunk fit onto the framing system. These effects could 
certainly have originated from shrinking of the shell around 
faying surfaces welded to the hull or possible a dynamic 
loading effect of the ocean waves. | 

The deformation of Figure25 (which is the authors 
version of the situation) illustrates that continuity of the 
structure between bulkheads must be considered in analysis. 
The extent of this significant structural section also requires 
a full consideration of shell curvature. One handicap of 
the gross section approach which is also obvious from Figure 
is that the localized stresses created by the superposed 
panel deformation would not be immediately determined. A 


final check of the solution should incorporate a method of 
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accounting for these localized stresses. 

The orthotropic rigidity coefficients (equation 1h) have 
been obtained and defined by the superposition of the susie 
laws of the isotropic shell and the mating gridwork frame. 

A better appreciation for these may be achieved by an examina- 
tion of these coefficients on an actual hull. Figure 26 
illustrates the magnitude and variation of the coefficients 
for the DD 931 class destroyer. The data for these plots 


is in Appendix D ° 
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constructed on DD 9 


Be Fig. 26 Orthotropic 





C. BOUNDARY CONDITIONS FOR USE wiTH 


PARTIAL DIFFERENTIAL EQUATIONS OF SHELL 


The solution of the developed partial differential 
equations (/5) will eventually require the stipulation 
of appropriate boundary conditions. These conditions will 
consist of a rational approximation of (a) the degrees of 
fixity existing at the bulkheads and at the junction of 
the curved shell and flat ships side and, (b) the magnitude 
and distribution of the loading which exists at the afore- 
mentioned boundaries. 
nui 3 

The degree of fixity which exists at the boundaries 
of the shell could indeed be the subject for a separate 
and absorbing analysis. Vedeler (1) suggests a method 
for accomplishing this. It would be of questionable value 
ovewer to pursue this subject in any depth in this thesis 
Since the immediate problem of determing a method of 
solution to the differential equation is ae major obstacle, 


For a first approximation then,. the following conditions 


are postulated. At the junction of the shell and the bulk- 
head, the fixity can be considered to be 100% clamped. The 
presumtion is a natural result when one considers that the 
structure is a ‘Shell continuous over supports (bulkheads) 
with continuous and equal loading on each side of these 


supports. An examination of the bulkhead N, stress 


ee 





distribution in Figure ]1] would also lead one to assume 
that a full clamping condition prevails. The fixity 
condition along the longitudinal boundary is a bit more 
elusive. Once more, as a first approximation, it will be 
assumed as being simply supported and capable of movement 
in the athwartship direction, The condition can be 


represented by Figure27. , 


a 





~\ 
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Figure 2 Degrees of Fixity at Shell Boundaries 
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Edge Loading; 

To begin the estimate of edge loading and-its dis- 
tribution, oe rel the bottom section and the total 
vertical force, due to the water (q) and the internal 





loading (w). 


where oe ye » No = a8 dese 


head for protection of freeboard 
plating. 


W = const. across section (*/r1+) 
dy = (44+ X z)cosh 

dy = %(ho+ z)cos¢ 

yi? = 2H cos¢ 

las = rLdd 

Fig. 28 Shell Loadin z = H/2 (1 + cos2h) 





total vertical force on hull Ve ) 
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An alternate expression for Vy could have been, 


. = fe - Wn - wt) AS, 


where Wy, Wy, are normal and tangential components of the 
internal loading w. dS is an increment of surface area 
projection on the horizontal plane, 

Substitution of w, = w cos0 

Wi = w sind 

dS, =  L cosf af 
and integration will give same results obtained previously 
for V+. 

Having now obtained the total upward force, a 
determination must now be made of the distribution of the 
reaction forces which support Vi and act on the shell 
boundaries. This is a statically indeterminate problem. 
In order to circumvent this obstacle and avoid an auxiliary 
analysis (whose accuracy would be greater than justified), 


the following approach is taken. 


| ~*&/, - la. te 





Fig. 29 Analo of Flat Plate Prismatic 
Loading to Hull Loading 
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Timoshenko (19 describes the behavior of a simply 
supported rectangular flat plate under a triangular 
prism distribution of loading, Figure 29. The portion: 
of this reference which is of interest is the nature of 
the reactive forces along the simply supported edges. 
Indeed it could be said that the vertical component of 
loading on the ships hull is a triangular prism distri- 
bution (with max value at the keel) if the hull were 
developed out into a flate plate, Figure 29. Because 
the ships loading is not on a flat plate but on a cylin- 
drical hull, no further resemblance or analogy will be 
made other than a consideration of the relative distribu- 
tion of the reactive forces and their relationship to the 
aspect ratio of the simple supports. Timoshenko tabulates 


V V 


y max » “x max for various aspect ratio of simple support 


boundaries. For the ship types being considered in this 
analysis, the pertinent ratio of B/L is approximately 1.3. 


This value gives 


Vy mex 2sf3 & aspect ratio 


Vx max 

or 
Sh 2s iyi 
VB 
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Therefore, ft 


eLV, t+ a BV, = Van 


2LV, + 2 BVg = HLY (“Tho + 3/41 H - kw ) 
x 


and, 
V ey (fame B 


VB 
Solving for the reaction forces gives, 


i, ae secs BHL® ¥ (ho , 307 a - ee 
2.7 L° +B2 h y | 


Vi = 0.5 BE (Tho 4 2TH by 
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Fig, 30 Edge Loading on Shell Boundaries 
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D. RIGIDITY COEFFICIENTS DATA 


The orthotropic rigidity coefficients of equation (1) 
can be written in the following form for y= 0.3. 


Extensional coefficients 
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TABLE 3 DD 931 RIGIDITY COEFFICIENTS DATA 
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DD 931 CALCULATED EXTENSLONAL RiGLLITY 
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COBFPFICIENTS AT BULKHEAD 
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TABLE 5 DD 931 RIGIDITY COEFFICIENTS (Actual) 
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TABLE 6 DE 1033 CALCULATED EXTENSIONAL RIGIDITY 


COBFFICIENTS AT BULKHEAD 
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TABLE 7 DE 1033 RIGIDITY COEFFICIENTS oe 
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